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Under heating, «-AlF3 undergoes a structural phase transition from rhombohedral to cubic at 
temperature T around 730 K. The density functional method is used to examine the T—0 energy 
surface in the structural parameter space, and finds the minimum in good agreement with the 
observed rhombohedral structure. The energy surface and electronic wave-functions at the minimum 
are then used to calculate properties including density of states, T-point phonon modes, and the 
dielectric function. The dipole formed at each fluorine ion in the low temperature phase is also 
calculated, and is used in a classical electrostatic picture to examine possible antiferroelectric aspects 
of this phase transition. 

PACS numbers: 64.70.Kb,71.20.-b,77.80.Bh 



I. INTRODUCTION 

The ionic insulator AIF3 has a number of known 
polymorphs 1,2 , which all convert irreversibly to the stable 
0J-AIF3 within the temperature range approximately 730 
to 920 K. Recent interest in this material arises due to its 
catalytic activity for dismutation and halogen exchange 
reactions^. 

Above its transition temperature (about 730 K)», a- 
AIF3 has the cubic perovskite structure AMX3, with the 
A cations absent (or the Re03 structure.) Aluminum 
plays the role of the M cation, and is surrounded by an 
octahedron of corner-shared fluorine atoms. At low tem- 
perature, the structure becomes rhombohedral, and this 
symmetry lowering can be characterized as a rotation of 
the fluorine octahedron about one of the three-fold axis 
of the perovskite cubic cell (or the a~a~a~ system in 
Glazer's tilt system^). As every adjacent octahedron ro- 
tates in the opposite sense, the wave- vector of this dis- 
tortion is (tt, it, 7r), and the unit cell becomes double the 
size of that in the cubic phase. 

We study the a-AlF3 structure by using pseudopoten- 
tials, a plane-wave basis, and the LDA method. After 
being tested, the pseudopotentials are applied in bulk a- 
AIF3. The total energy surface is examined in the struc- 
tural parameter space where the cubic phase is compared 
with the rhombohedral phase found at the deepest nearby 
minimum of the energy surface. 



II. TESTING DIFFERENT 
PSEUDOPOTENTIALS 

The Troullier and Martins method^ is used to generate 
the pseudopotentials. Since there are different possibil- 
ities of choosing cut-off radii r c /, as well as how many 
valence orbitals to include in the pseudopotentials, we 
test the pseudopotentials by doing LDA self-consistent 



iterations for a crystal or molecule prototype and com- 
pare the results with FP-LAPW— all-electron calculations 
(full potential linearized augmented plane- wave 8 ). This 
provides guidance on the choice of energy cut-off of the 
plane-wave basis, and also the choice of local potential, 
in the a-AlF3 calculation that follows. The results of 
these tests are shown in Table [I] and Table [HI for alu- 
minum and fluorine separately. For the Ceperley-Alder 
exchange-correlation energy^, the Vosko-Wilk-Nusair 10 
parametrization is used. 

The test of the aluminum pseudopotential is performed 
in bulk Al metal, not only the observed fee structure, 
but also the sc and bee structures, where comparison can 
still be made between the all-electron method and LDA. 
The results show that the pseudopotential for the Al 3+ 
ion should be capable of describing the unoccupied 3d 
orbital. Otherwise, even the large r c i (I = 0, 1 for s and 
p orbitals) needed to fix the lattice constant in the fee 
aluminum are insufficient to remove the r c ;-dependence of 
the lattice constant and bulk modulus. Due to the more 
attractive nature of the d pseudopotential, including the 
3d orbital reduces the dependence of these two quantities 
on r c i and the choice of local potential, and improves the 
agreement with the all-electron results. Better agreement 
is also seen in comparisons for sc and bee structures of 
Al metal. 

The fluorine pseudopotential is tested for the F2 
molecule. Unlike the low plane-wave cut-off needed for 
aluminum (less than 40 Rydberg), more than 90 Ryd- 
berg is needed in the F2 molecule. For bulk Al metal we 
use a k-point mesh of 12 x 12 x 12 that contains 56 special 
k-points, while for F2 we do a one k-point calculation. 

Results in the following sections are calculated using 
pseudopotentials (Al^, F^) as shown in the captions of 
the Tables. This set of pseudopotentials gives the lattice 
constant (a — 5.02A), and bulk modulus (B = 146.6 
GPa) of cubic a-AlF3 in good agreement with all-electron 
results (a = 5.043A, B = 150 GPa.) 
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TABLE I: Aluminum pseudopotential tests on bulk Al. The 
first column shows the cut-off radii in units of Bohr radius 
(do) for the valence orbitals (1=0,1,2). Tests with r C 2 omitted 
do not include the d orbital. The following columns show the 
use of core correction, the choice of local potential, lattice 
constant, bulk modulus and the total energy per atom with 
respect to the fee structure. Values are shown in atomic units 
(e = electron charge, ao = Bohr radius.) AE is defined to be 
zero for the fee structure, while the total energy differences 
between fee and other structures can be compared with all- 
electron results. The pseudopotential with a ' is chosen for 
the calculations of a- AIF3 . 
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III. RESULTS FROM LDA 

A. k-point Test and Density of States 

A k-point mesh of 3x3x3 is used for the Brillouin zone 
summation for a-AlF3. This is tested for cubic CV-AIF3, 
where a finer mesh of 4x4x4 gives a barely different en- 
ergy vs plane wave cut-off curve, with energy uncertainty 
less than 1 meV. This mesh gives six special k-points after 
symmetrization (according to the point group of the 



TABLE II: Fluorine pseudopotential tests on F2 molecule. 
The last column indicates the plane-wave cut-off needed to 
obtain the convergent result. The pseudopotential with a ' is 
chosen for the calculations of Q-AIF3. 



r c o/r c i 


core 


local 


bond 


force 


cut-off 


(a ) 


corr. 




(a ) 


const. (e 2 /ao) 


(Ryd.) 


1.4/1.6 


110 


P 


2.780 


0.370 


90 


1.4/1.45 


yes 


P 


2.710 


0.408 


90 


1.3/1.3 


yes 


s 


2.599 


0.384 


over 100 


1.3/1.3 t 


110 


s 


2.603 


0.376 


90 


1.2/1.2 


yes 


s 


2.617 


0.387 


over 100 


all-electron 






2.632 


0.376 




experiment 






2.669 


0.302 





crystal), over which we include at least 36 states in the 
truncated Hamiltonian diagonalization processes. Since 
in the rhombohedral unit cell of CV-AIF3 we consider the 
48 valence electrons from aluminum 3s 2 3p 1 and fluorine 
2s 2 2p 5 , these electrons will occupy the 24 lowest lying 
bands. A test including more states, 48 states for exam- 
ple, shows that in the process of iterative diagonalization, 
using 36 states gives convergent answers for the eigenval- 
ues of the lowest 24 bands. Thus we include 36 states in 
the calculation of the minimum on the energy surface. To 
make the density of states plot, we include 50 states to 
better describe the empty states from aluminum orbitals. 
As shown in Fig. ^ LDA gives an insulator band-gap of 
8 eV. We choose 90 Rydberg to be the energy cut-off for 
the plane-wave basis. This choice gives the same energy 
difference between the two structures as that from the 
choice of 140 Rydberg (with error within 1^2 meV. 
The two structures are the cubic and the rhombohedral 
structure at the minimum.) 




Energy(eV) 



FIG. 1: Density of states plot of Q-AIF3 in the rhombohedral 
ground state, using pseudopotential (A\^ , F^). The rhombo- 
hedral ground state has unit cell edge 4.76A, 6=57.5° and 
5=0.098. This calculation includes 50 states and uses a k- 
point mesh 6x6x6. The temperature broadening is 54 meV. 
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B. Crystal Structure 

We use LDA to search for the ground state in the rhom- 
bohedral symmetry class i?3C seen experimentally in a- 
AlF^. There are three free parameters: (1) the lattice 
constant a, (2) the angle 9 between any two of the rhom- 
bohcdral lattice vectors, and (3) the inner parameter 5, 
which describes the bending of the Al-F-Al bond. To de- 
scribe these three parameters, we first start from the cu- 
bic structure with lattice vectors (26,0,0), (0,26,0) and 
(0,0,26), where 6 is the Al-F bond length. A rhombo- 
hedral unit cell is formed by picking new lattice vectors 
ai = (0,26,26), a 2 = (26,0,26), and a 3 = (26,26,0). 
The angle 9 defined as 9 = cos _1 (ai ■ «2 / 1 a 1 1 1 ^2 1 ) is 60° 
in this case. This new cell is then elongated along the 
[111] direction. Consequently, the strains €23 = £31 = £12 
become nonzero, and 9 7^ 60°. Viewed from the [111] 
direction, in the aluminum plane the ions form a hexago- 
nal pattern, and the nearest Al-Al distance on this (111) 
plane is the lattice constant a. This rhombohedral struc- 
ture can also be described in a hexagonal setting, with 
the following c/a ratio 



c 
a 




1 + 2 cos 6 
1 — cos 9 



(1) 



The cell shape and cell volume V are determined by a 
and c/a: 



V 



1 + 2 cos ( 



2\/2V 1- cos 9 2V^' 



(2) 



while the relation between the lattice constant a and the 
length of the rhombohedral lattice vectors is: 



a/\3i\ = v/2(l - cos6>) 



(3) 



There are two aluminum ions and six fluorine ions in each 
rhombohedral unit cell. The fluorine ions sit at the 6e 
sites ((x,x + j), (x,x + \, |), etc.), and (5=x-0.75 is 
the deviation from x—0.75 where the Al-F-Al bond angle 
is strictly 180°. The distortion caused by S alone defines 
the octahedron rotation angle u> about the [111] axis: 



tanw = 2v / 3<5, 



(4) 



On the other hand, the c/a ratio and the lo have the 
following relation: 



V6a 



i + C 



coscj 



(5) 



where £ describes the flattening (£ < 0) and elongation 
(C > 0) of the octahedron along [111]. 

We first assume cubic symmetry by fixing 9 = 60° and 
6 = 0. LDA gives the lattice constant a = a* = 5.Q2A 
for the cubic phase. With the lattice constant kept fixed 
at the value a = a*, 6 is then relaxed to give minimum 
energy at (a, 9,5) = (a*, 60°, ±0.04). As shown in Fig. 



13 the relaxation of S finds minima at £=±0.04 and a 
maximum at 6 = 0, with an energy difference of 14.4 
meV, showing that in our T=0 energy surface study, the 
unbent Al-F-Al bond, and therefore the strictly cubic 
structure, is not a metastable solution. 

The three parameters are then relaxed in turn to ap- 
proach the deepest nearby minimum. In fact, in the pa- 
rameter space, there is a path that keeps the fluorine 
octahedra rigid and leads to the region near the over- 
all minimum. A Rigid octahedron means that the Al-F 
bonds are fixed at the bond length found by LDA in the 
cubic phase 



(0) 



2V2 2^2 cost 



and that the angles between neighboring Al-F bonds are 
always 90 degrees (( = 0). With these two restrictions, 
only one parameter S remains and defines the rigid octa- 
hedron path. The energy along this path is also shown 
in Fig. 12 The overall minimum locates not far from 
the minimum on this path. It has an Al-F bond length 
only 0.3% longer, a pretty small octahedron elongation 
factor ( = 0.25%, and the total energy 93 meV lower 
than the cubic phase. Table IIIII shows the comparisons 
between LDA and experimental data 4 . While no T = 
experimental numbers are available, one can only make 
reasonable comparisons for the deviation of bond length 
and £. The rhombohdral phase observed in the exper- 
iment has an Al-F bond length 0.43% longer than the 
value observed in the cubic phase, and ( = —0.102% 
(meaning the octahedron is slightly flattened.) This in- 
dicates that both the LDA and the experimentally found 
rhombohedral phases are close to the rigid octahedron 
path. 



TABLE III: The structure parameters found by LDA are com- 
pared with experimental values- here for different phases. The 
total energy (per two formulas of AIF3) difference between the 
two phases is also shown here, where the total energy of the 
cubic phase is set to zero. The values for the rhombohedral 
phase are calculated at the minimum of the energy surface. 
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Phonon Modes 



Previous workers^ investigated the phonon spectrum 
of the cubic phase using the generalized Gordon-Kim 
method 12 . Density functional theory was applied to con- 
struct the charge distribution and polarizability of ion. 
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FIG. 2: The relaxation of 8 from the cubic structure is shown 
in plot (a). The curve of total energy (per two formulas of 
A1F 3 ) vs 5 at 9 = 60.0° and a = 5.02A shows that the cubic 
phase is sitting at an energy maximum. In (b), the curve 
follows the rigid-octahedron path, and leads to the region near 
the overall minimum in the parameter space. The dotted- 
dashed line shows the energy of the overall minimum. 



An approximate crystal energy is then calculated from 
the ion results, and is used to examine the crystal lat- 
tice dynamics. However they did not find an instability 
of the cubic phase, which is contradictory to the experi- 
ment. Here we use full LDA to examine the energy sur- 
face around the rhombohedral minimum, and extract the 
soft phonon modes directly. The rhombohedral minimum 
described in the previous section, and the dipole aspect 
discussed in Sec. IIVI both show the instability of the 
cubic phase. 

At the r point of the rhombohedral phase, the ir- 
reducible representations are those of the point group 
D^d (note that the space group of the rhombohedral 
phase is i?3C, which is nonsymmorphic.) Among the 
12 r point modes provided by the fluorines, namely, 
A lg 9 2A 2g 8 3E g A lu ® 2A U e 3E U , the four modes 
Ai g ffi 3E g are Raman active. Viewed from the cubic 
phase, the octahedron rotation that governs the struc- 
tural transition is the R 5 mode (in Kovalev labelling) at 
the zone boundary (tt, it, 7r)lA. After distortion, the R5 
mode gives rise to the zone-center Ai g mode and one 
of the three E g modes. Raman experimentsi^i^ showed 
that the Ai g and one E g are the soft phonon modes below 
transition. The Ai g is the Al-F-Al bond angle bending 
mode. Using the fluorine atomic mass, and the energy vs 
5 curve, the Ai g frequency is calculated to be 205 cm -1 . 
We also calculate the other restoring coefficients of the 
three E s modes, and the two Raman-inactive A2 g modes. 
Their frequencies are listed in Table. IIVI 



TABLE IV: Some of the F point phonon modes frequencies 
calculated by LDA. direct comparison with experiments are 
not available since previous observations are done at RT and 
above. However, in the experiment it is found that the two 
high energy E g modes only weakly depend on temperature. A 
crude extrapolation using uj 2 — A(T — T c ) (from Curie-Weiss 
law and Lyddane-Sachs- Teller relation) and the RT experi- 
mental data™ gives the frequency values at T = 0. 
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D. Dielectric Function 

The electronic part of the dielectric tensor e(q — 0, ui) 
is also calculated. e xx — e yy 7^ e zz is expected since we 
have chosen lattice vectors such that x and y are perpen- 
dicular to the three- fold axis. The imaginary part of the 
dielectric function can be directly obtained by calculating 
all direct inter-band transitions 15 : 



e a (w) = (— ) 2 ^^|(Vn )fc b1^ )fe )| 2 /»(l-/ m ) 
muj * — ' * — ' 

k m,n 

x 6{E m - E n - Tllu), (7) 



where p is the momentum operator, and k runs through 
all k-points in the Brillouin zone allowed in a unit vol- 
ume. The band indices are m and n, and f m ,k are 
the occupation number of the n th and m th states at the 
k th k-point. The real part of dielectric function is ob- 
tained using the Kramers-Kronig relation. The result at 
T=0 is given in Fig. [21 but we do not know of any exper- 
iment to compare with. However this calculation does 
neglect the non-local potential effecfei^ and the local field 
effect, and follows the pseudo-wave-functions that are 
smoothened in the ion core regions. 



IV. INNER PARAMETER S AND DIPOLE 
FORMATION 

The formation of a dipole plays an important role in 
the energy difference between the two a-AlF3 structures. 
When the inner parameter 5 is non-zero, each fluorine 
atom develops a dipole (Aluminum atoms sit on inversion 
centers and therefore cannot have dipoles). Starting from 
the cubic structure, the distortion from 5 alone defines 
the octahedron rotation angle u> about the [111] axis (see 
Eq. UJ and already lowers the symmetry to R3C. The 
parameter 5 is irrelevant to the c/a ratio, but it does 
flatten the octahedron (see Eq. [SJ). It also elongates the 
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FIG. 3: Plots of dielectric functions. This calculation includes 
50 states and uses a k-point mesh 6x6x6. 



Al-F bond length, and introduces a polarization to each 
fluorine. The energy gain from the fluorine dipole-dipolc 
interaction drives the structure off-cubic. This is shown 
in FigEl where the cubic structure is found to be not a 
metastable solution. 

A naive picture of ionic solids has spherical electron 
charge clouds of total charge Qi around the i th ion. In 
reality, charge clouds are distorted. For example, in high 
T cubic a-AlF3, the F ions are noticeably prolate when 
examined in a (lOO)-plane charge contour calculation^. 
This distortion is still evident in the low T rhombohedral 
phase, as shown in Fig. 0] The value of the charge Q^, 
on the other hand, is not uniquely definable. By creating 
a sphere centered at each fluorine atom, with a radius 
that only allows neighboring spheres in the rhombohe- 
dral phase to touch each other at one point, we define 
a volume to examine the charge and the dipole moment 
of the fluorine. In the cubic phase, such a sphere con- 
tains charge -0.675|e|, while it contains -0.68|e| in the 
rhombohedral phase (with a slightly different radius.) 
Using the same sphere, we compute the local induced 
dipole moment p ind of the fluorine ion (which also has 
no unique definition 19 .) As expected from symmetry, 
find = in the cubic phase, whereas p; n d = 0.103|e|A in 
the rhombohedral phase, pointing from the fluorine posi- 
tion of the unbent Al-F-Al bond towards the actual dis- 
torted fluorine position. Note that the displacive dipole is 
Pdis = 0.315|e|A, defined by the fluorine charge -0.675|e| 
and the displacement 0.467A of the fluorine ion when the 
Al-F-Al bonds bend. By symmetry, the other fluorine sit- 
ting at the opposite side of the octahedron has its dipole 
pointing exactly in the opposite direction. Therefore the 
material could be called antiferroelectrio^S. 

The dipoles and charges can as well be determined by 
doing integration in one of the equivalent polyhedra that 
surround the F~ ions and partition the whole space (the 




FIG. 4: Charge density plot of the rhombohedral phase on 
a (100) plane. The [100] direction is ill-defined since it is 
inherited from the cubic phase, but the plane shown here 
does contain the aluminum ions (labelled in solid dots) that 
form the (100) plane in the cubic phase. The fluorine ions are 
slightly off the plane shown here, and display clearly a prolate 
charge distribution. 



Al 3+ 's sit on shared vertices of the F polyhedra.) In 
this way, the ions are charged as Al 3+ and F~, while 
p- md = 0.0963|e|l and p d[s = 0.467|e|A. 

It is worthwhile to see if a simplified picture with these 
numbers helps to understand the energy. In LDA, the 
energy includes the exact electrostatic or Hartree energy 
from the electron charge clouds, as well as other quan- 
tum effects (energy of derealization, exchange and cor- 
relation.) Purely classical electrostatic models (plus hard 
core repulsion), on the other hand, provide a simple but 
useful view. We shall see how much classical electrostatic 
energy comes from the charge and if this helps to explain 
the stability of the structural deformation. 

In a simplified picture where only electrostatic energy 
and dipole formation energy (Ei in Eq. |H1 where a = 
0.858A 3 — ) are considered, the total energy is 
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This form leaves out the quantum effects and the effect 
from covalent bond angle. We find that E%, which is 
purely ionic electrostatic (including displacive dipoles), 
is less negative in the rhombohedral structure than in the 
cubic structure (i.e., if the same ion charge assignment is 
considered in both cases.) However, when the induced- 
dipole-related terms {E\ + E 3 + E 4 ) are also considered, 
the energy is lowered. By assuming point charges and 
point dipoles, we use the results from the polyhedron 
method to illustrate this picture: E 2 per two A1F 3 for- 
mulas in the rhombohedral phase is 1 eV higher than that 
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in the cubic phase; but E% and E4 bring the energy E 
down, with the cost of Ex , to 0.97 eV lower than the cubic 
phase (which has only E 2 ) This suggests that the effect 
of induced dipoles is more than enough to compensate 
the energy loss from the structural distortion. 

V. CONCLUSION 

We report an LDA study of bulk a-AlF 3 . By examin- 
ing the T=0 energy surface for structures of the phases 
on both sides of the transition, we find the structural 
parameters to agree with previous experiments, and the 
cubic phase not to be a metastable solution. Using the 
result of LDA, the density of states plot and the dielec- 
tric function are provided. At the T-point, the predicted 
Ai g soft phonon mode and E g modes are compared with 
previous Raman experiment, while two other A 2g modes 
are predicted. We look at the charge and dipole moment 



at each fluorine ion, and use these quantities to calculate 
the classical electrostatic energy. In the antiferroelectric 
distortion which accompanies the structural transition, a 
classical calculation shows that the electrostatic energy 
gain from dipoles is more than enough to compensates 
the energy loss from the ion-array deformation. 
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